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We present a theory for tunneling spectroscopy in a break-junction semiconductor device for ma- 
terials in which the electronic conduction mechanism is hopping transport. Starting from the con- 
ventional expression for the hopping current we develop an expression for the break-junction tunnel 
current for the case in which the tunnel resistance is much larger than the effective single-hop re- 
sistances. We argue that percolation like methods are inadequate for this case and discuss in detail 
the interplay of the relevant scales that control the possibility to extract spectroscopic information 
from the characteristic of the device. 



72.20.i, 73.40.Gk, 73.40 Lq 



I. INTRODUCTION 



Tunneling spectroscopy is one of the most widely used 
tools for obtaining information on the electronic struc- 
ture of solids^'^. If the tunneling takes place between 
electrons of the same energy (elastic or resonant tunnel- 
ing) the tunneling current is a convolution of the densities 
of states (DOS) of the contact materials^. Thus, the DOS 
of one material can be extracted if the DOS of the other 
material is known. In the case of inelastic tunneling the 
electronic transitions occur with emission or absorption 
of phonons. In this case the tunneling current becomes 
also sensitive to the phonon DOS'*. 

Applying tunneling spectroscopy to doped semicon- 
ductors in the hopping regime at low temperatures has 
proven to be very useful in revealing the influence of the 
electron-electron interaction on the DOS of the impurity 
band in the meV range. Such experiments have been per- 
formed with conventional metal-barrier-semiconductor 
contacts^'^ as well as with mechanically controllable 
break junctions^. Both methods reveal the Coulomb gap 
in an impressive way. 

The theoretical interpretation of break-junction ex- 
periments in which both contacts consist of a doped 
semiconductor^ is more involved than that of a metal- 
insulator-semiconductor structure^. As the physical na- 
ture of the inelastic tunneling transitions between the 
contacts appears to be similar to those between the lo- 
calized states in hopping transport the question could 
arise whether all transitions in question are a part of a 
global disordered resistor network which would have to 
be treated by conventional percolation methods. The as- 
sumption that this is indeed the case has, for example, 
been the basis of the arguments in Ref. [ 9] . 

At present it does not seem clear, whether in the hop- 



ping regime the traditional method for the calculation of 
the tunneling current due to Bardeen^ is applicable or 
has to be replaced by a more sophisticated version. 

Even the question, whether or not information on the 
electronic DOS can be extracted from break junction tun- 
neling experiments is discussed controversial in the liter- 
ature. In Ref. [ 9] it is claimed that the tunneling current 
should exhibit strong mesoscopic fluctuations as a func- 
tion of the applied voltage. Indeed, if this would be true, 
little information on the global electronic DOS of the ma- 
terial could be obtained. On the other hand, other the- 
oretical arguments"'^^ lead to the conclusion that, under 
certain conditions, information on the electronic DOS can 
be obtained. Furthermore, neither in experiments^ nor 
in simulations of tunneling between materials with local- 
ized states^* strong fluctuations of the tunneling current 
as a function of the voltage are observed. 

Therefore we address this issue here from a fundamen- 
tal point of view. We start from the rate equations for 
hopping transport and derive an expression for the tun- 
neling current on the basis of these equations. The fur- 
ther development of the theory exploits the fact that the 
typical transition probabilities in hopping transport are 
orders of magnitudes larger than those across the junc- 
tion. This is due to the fact that the junction separation 
is large compared to the characteristical hopping length 
in the bulk and due to the fact that the inverse of the vac- 
uum tunnel parameter k is much smaller than the local- 
ization length a~^, which is the tunnel parameter of the 
hopping transitions. Therefore a separate equilibrium is 
established in both contacts with two separate chemi- 
cal potentials the difference of which is controlled by the 
bias voltage. The resulting expression for the current is 
similar to the conventional expression for tunneling spec- 
troscopy. Using this expression we discuss under which 
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conditions concerning the length scales and the electron- 
phonon coupling simplifications of the current formula 
occur. 



II. THE TUNNELING CURRENT 

In our derivation of an expression for the tunnel current 
in a break junction made of a lightly doped semiconduc- 
tor we first recapitulate the general theory of hopping 
transport in an impurity band of a bulk semiconductor, 
then we consider the case of two parts of a break junc- 
tion separated very far from each other. These steps will 
then enable us to study the break junction under realistic 
conditions. 

In the standard theory of hopping transport^^'^"^ the 
impurities are assumed to provide localized electronic 
states at sites Rm,R„ with localization length and 
characteristic energies Charge carriers (electrons 

or holes, dependent on the type of doping) perform 
phonon-assisted tunneling transitions (hopping transi- 
tions) between these states under the influence of an 
externally applied electric field E. The interaction be- 
tween the charge carriers is assumed to be Coulomb-like 
(Hubbard-interaction effects are ignored). For definite- 
ness we assume the dopant to be n type, i. e. we con- 
sider electrons with charge q = —e {e = \q\ is the ele- 
mentary charge). If the interaction is treated in Hartree- 
Fock (HF) approximation (i. e. if many-particle jumps 
are ignored) the dynamics of the charge carriers can be 
described by the simple rate equation 



dt 



^{l - Pn)W^n] (1) 



Here is the probability to find a charge carrier at site 

m, 

(3 

Wn,n = l^i\Vnm\) exp{-2aRnm + "(Km " \Vnm\)} (2) 

is the transition probability for a hop from the site n 
to the site m, (3 = l/ksT, Rnm = |Rram| is the distance 
between the sites (R„m = R„-Rm), and Vnm = V„-Vm, 
where 



(3) 



is the energy of a charge carrier at site m, (e is the di- 
electric constant of the host). z^(|Kjm|) is the spectral 
function which describes the influence of the electron- 
phonon coupling ("attempt-to-escape frequency"). The 
quantity i'(lKmil) characterizes the ability of the phonon 
to induce the transition. Since phonons with different en- 
ergy can interact with localized electrons differently well 
this frequency is a function of the energy transferred in 
one hop. 



As a model for a break junction we now consider two 
samples made of a lightly doped semiconductor which are 
separated by a distance / (sec Fig. 1). Wc distinguish be- 
tween sites situated on the left and right with respect to 
the separation ("left sites" and "right sites") and denote 
the labels of the left sites with lower case letters m,n 
and those of the right sites with upper case ones M, A''. 
Since the junction is a break junction both samples have 
the same densities of states in the absence of the electric 
field, that is Nl{V)\ E=0 = Nii{V)\E=0- 

If both samples are separated very far from each other, 
there are no transitions between left sites and right sites. 
Therefore, Eq.(l) is valid for each of the samples sepa- 
rately. Since both samples are in equilibrium, the solu- 
tions to the transport equations are simply Fermi func- 
tions with two different chemical potentials and ^n. 
For the left sample we have, for example 



Pn 



fm — 



1 



l + exp{/3(K.-ML)}' 



(4) 



If wc now decrease the sample separation, tunneling 
transitions between the left and the right sample become 
possible with transition probabilities 

f P 

WmM = t^{\VmM\) exp < -2KlmM + ^(HraM - |KnM|) 

(5) 

The rate equations for the occupation probabilities 
Pm{t) and pM(t) acquire additional terms of the form 

Pm(1 — Pm)WMm- It is important to note that the 
tunnel transitions betwc;en the left and the right sam- 
ple do not contain the localization length a but the vac- 
uum tunneling parameter k between the two bulk sam- 
ples. Kr^ is of the order of a few Angstroms, whereas 

typically takes values around 10 nm. Moreover, the 
site separation RmM is replaced by the effective tunnel- 
ing length IrnM = I + SlrnM ■ ^UnM if' ^ Small corrcction to 
the tunneling distance between the two samples taking 
into account the different wave function amplitudes for a 
given pair (m, M) of localized states. 

If one knows the solution to the transport equations 
and a proper method to deal with the configuration aver- 
age one can calculate the current by averaging the quan- 
tity 



J = 



"dT ^ 



M 



At 



(6) 



Here f2 is the total volume which contains the sites con- 
tained in the summation in Eq.(6). 

Wc now distinguish between two fundamentally dif- 
ferent situations which can be controlled by the experi- 
mentally adjustable parameter I. If Z is small enough, a 
common equilibrium between the two samples can be es- 
tablished. In this case the standard methods for solving 
the transport problem in a disordered hopping system. 
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namely Mott's optimization-'^^, percolation theory-'^^-'^ or 
the effective-Medium approximation^^ can be applied. 
However, if 2kI becomes appreciably larger than the ex- 
ponents of the dominant bulk transition probabilities (2) 
(which is already the case if I becomes larger than a few 
nm) , a separate equilibrium is established in the two sam- 
ples, before a tunneling transition can take place. In this 
situation, which is the one we are interested in, the tun- 
nel transitions are not a part of a percolating network 
because the charge carriers that cross the break junction 
cannot optimize their paths. In fact, the tunneling ex- 
periments reported in Ref. [ 7] arc performed in such a 
way that the resistance of the tunneling contact is by 
orders of magnitudes larger than the equilibrium bulk re- 
sistance Req, SO that the current is determined by jumps 
across the contact. Furthermore, care was taken to adjust 
the distance I in such a way to make sure that the resis- 
tance Rtu across the junction did not show the temper- 
ature dependence of the bulk resistance Req oc e^^°/^^^ , 
X « 0.5 which would indicate that a common thermal 
equilibrium of the two contacts would exist. In the range 
Rtu ~ 10^ Req this regime was reached. It is this regime 
in which the Coulomb gap was observed. Accordingly, 
a charge carrier can hardly optimize its path by return- 
ing across the junction. We therefore make use of the 
separation of time scales in the present situation. 

The time for a tunneling transition can be estimated 



as 



t 



tU ' 



exp{2K/} 



whereas the equilibration time is roughly 
exp{(To/T)-} 



eq 



(7) 



(8) 



so that we have Rtu /Req = Uu/teq- As mentioned above 
this ratio is of the order of 10^ in the tunneling experi- 
ments, so that one can be sure that the two sample equili- 
brate separately before the tunnel transitions take place. 
In this situation the quantities Pm and Pm can be con- 
sidered as Fermi functions as before, and we obtain for 
the tunneling current 



j = -^(1 - e-'^(''«-^-)) R„M fm{l - fM)W, 



mM- 



%M 



(9) 



Let us now discuss how to perform the configuration 
average over the current (9). First we note that tun- 
nel transitions which lead "upwards" in energy are 
extremely unlikely. So we are left with WmM = 
J^(lKraMl) cxp{— 2KZmM}- Due to the exponential depen- 
dence small fluctuations of ImM lead to large ones of 
WmM- Despite this fact the sum in (9) itself is not 
strongly fluctuating, since it is a sum of many terms 
which are nearly statistically independent of each other. 
The sum has the structure (j = jxSx) 



jx — ^ ^ jmM ■ 



(10) 



The currents jmM and jm'M' in Eq.(lO) are only corre- 
lated if either m = m' or M = M', or both equalities 
hold. Therefore, deviations from the configuration aver- 
age satisfy the relationship 



(01)) 



oc 



1 

N' 



(11) 



Here ((jD) = (j1) — (jx)^ is the standard deviation, (. . .) 
symbolizes the configuration average, and N is the num- 
ber of initial sites, which was assumed to be same as 
the number of final sites. The argument proves that the 
configuration average can be used for sampk;s with suf- 
ficiently large contact area. The situation for the exper- 
iments of Ref. [ 7] is depicted in Fig.l. In this experi- 
ment the contact area is of macroscopic size (Ly ?s0.8 
mm,/^ «0.2 - 0.4 mm). Accordingly, N is a large number 
(Although the total number of initial sites in the contact 
area is of the order of 10^ the quantity N must he some- 
what smaller, since only pairs in a strip of width eU + kT 
in energy space contribute to the sum. The actual num- 
ber depends on the width of the impurity band, which 
is not exactly known). Accordingly, we expect that jx is 
also a self-averaging quantity for this experiment. This 
assumption is confirmed further by the fact that the fluc- 
tuations of the tunneling current in the experiment were 
smaller than 10 percent, whereas huge fluctuations would 
be expected outside the range of applicability of the av- 
eraging procedure. 

For calculating the configuration average we use the 
densities of states Nj^jjiiy, R). In the presence of an elec- 
tric field they are quantities which depend on the energy 
V and on the position vector R. The latter dependence 
describes the space charge region, that is the region, in 
which the charge carrier density in the presence of the 
field differs from that in the absence of the field. This 
dependence is only negligible deep in the interior of the 
sample. There wc have N^/n = Ni/ij{V) independent 
of R. Performing the configuration average by means of 
the densities of states we obtain the expression 

(j) = -^(l-e-^(''™)) 

X j dR j dR' j dVdV'{R-R')NL{V,R) (12) 

X Nr{V', R')fR{V'){l - fUV))Wi\R - R'\,V', V) 

for the configuration averaged current. 

We would like to emphasize that the same averaging 
procedure can also be applied to the conductivity of the 
bulk^^. However, in this case a different expression for 
the current has to be used, which takes into account that 
the particle optimizes its path through the sample. Doing 
so, it often returns to its initial site. Therefore, the distri- 
bution functions become functions of the transition prob- 
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abilities, so that /,„ Pm{{Wrnn})- The latter quan- 
tities are calculated from the Miller-Abraham random 
resistor network. Percolative aspects of hopping trans- 
port are included if the statistical correlation between 
the transition probabilities in the distribution functions 
is taken into account, as it is done in the effective-medium 
theories. In a tunneling experiment, however, the situa- 
tion is different. A particle, that has managed to cross 
the junction, never returns to its initial site to look for 
a better path, since every hop it can perform in its new 
surrounding is easier to perform than a hop across the 
junction. Accordingly, the particles equilibrate in their 
new surrounding on a time scale which is small com- 
pared to the tunneling time. Mathematically, this fact is 
expressed in that the occupation numbers are indepen- 
dent of the transition probabilities, so that the quantities 
jnm in Eq.(lO) depend only on n and on m, but not on 
any other site. The latter fact is a consequence of the 
separation of time scales. 

Eq.(12) is all what the kinetics tells us. If we want to 
simplify this equation further we have to use additional 
knowledge on the samples, that is, about the length scales 
present in the system. Such scales are the localization 
length a^^ , the tunneling distance I, the screening length 
le, and the sample length Lx- Different relationships be- 
tween these length scales yield different expressions for 
the current, as discussed further below. 



III. METAL-LIKE CONDITIONS 

The expression (12) takes a particular simple form in 
the case of metal- like conditions, i. e. le I, ^ 
le, and le -C a^^. The first condition means that we 
can use the approximation — = e.El = eU. Here 
U is the voltage applied to the sample. If the second 
condition holds there is a region in the samples in which 
the densities of states are independent of R. If the third 
condition applies mainly sites outside the space charge 
region contribute to the integral (12). so that the DOS 
measured is the densities of states in the bulk. 

In the limit of strong disorder it is unlikely to find 
neighboring sites on opposite sides of the break junction 
with the same energy. Therefore, mainly inelastic tran- 
sitions are relevant in this limit. Furthermore, as men- 
tioned before, jumps upwards in energy can be ignored 
at low temperatures since there are plenty of accessible 
sites, which can be reached by jumps down in energy 
space. Accordingly, we obtain^'* for Pell > 1 



— I dR„ / dRM(Rm - R.M)exp{-2K;ZTOM} = -H^e^;. 



(ix) =ewj dVdV Nl{V)Nr{V') 

xe{v' - v)e{i^R - v')e{v - fiLM\v - v'\), (13) 



where we have defined 



(14) 



Here 17l {^r.) is the volume of the left (right) sample, 
over which the integration takes place. Since the transi- 
tion probabilities are exponentially small quantities with 
respect to the tunneling length k^^ the range of integra- 
tion penetrates only over a distance of the order of a few 
times of into the sample. Thus the relevant volume 

is of the order of Al, where A = Lyl^ is the area of the 
cross section of the break-junction. 

Since the junction is a break-junction the densities of 
states Nl and A''^ agree with each other if the electric 
field is switched off, as noted above. Therefore, they 
differ only in the position of the zero point of the energy 
axis if the electric field is switched on. That is Nr{V) = 
Nl{V — eU). Accordingly, we obtain 

IXL+eU V' 

Ux) = eW J dV J dVNL{V)NL{V' - eU)u{V' - V). 

(15) 



IJ-L HL 



If it were not for the function uiy' — V), which describes 

the energy dependence of the electron-phonon coupling, 
we would now have a tool for extracting information on 
the density of localized states. If this energy dependence 
is not known one might have difficulties in interpreting 
inelastic tunneling spectra^°. 

If we assume the deformation potential approximation 
to hold and that a Debye model for the phonons describes 
the situation adequately well, the energy dependence of 
V is known and the integrals in (15) can be evaluated. 

In deformation potential approximation the function 
v(E) takes the form 



v{E) = vq 



\E\ 



(16) 



where s is the velocity of sound, and uq is a constant^^. 
This approximation takes into account that the overlap 
between the phonon-wavefunction and the wavefunctions 
for localized electrons decreases rapidly, if the phonon- 
wavelength becomes small compared to the localization 
length. 

Of particular interest is the situation in which the DOS 

shows a pseudogap centered at the Fermi energy, as it 
is the case in the presence of a Coulomb gap at finite 
temperature. In this case the DOS has the structure 



Nl{V) = No + n^\v - ^XL\ 



(17) 



where 7^2 for three-dimensional systems at zero tem- 
perature and A^o vanishes at zero temperature^^' Using 
(16) and (17) in (15) we obtain (see appendix) 
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Ox) <x \U\ 



3+C 



(18) 



for e\U\ <^ hsa, where C ^ depends on the parame- 
ters appearing in (16) and (17). Therefore, the data for 
the tunnehng conductance appear to scale to zero in an 
experiment, which is performed in the regime f3eU > 1. 
Since the apphed voltages are very small in the regime 
(3eU <C 1 we expect that this is also the behavior, which 
would be observed in experiments. However, wc would 
like to stress that the true value of the tunneling conduc- 
tance at zero bias is non-zero. To calculate the derivative 
of the current at zero bias we use Eq.(12). Doing so, we 
obtain 



d(jx) 
dU 



u=o = e^W{kTf[N^Joa + 2NoN^{kTyJo^ 



where 



J, 



7A 



4/ 



\x - y\ exp(-|a: - y\/2) 
cosh(a;/2) cosh(y/2) 



(19) 



(20) 



In the appendix also results for e\U\ ^ hsa arc presented. 

According to the Eqs.(18) and (19) the tunneling con- 
ductance scales to zero with decreasing U for f3eU > 1 
and approaches a constant at f3eU <C 1. The zero bias 
tunneling conductance itself increases at least quadrati- 
cally with increasing temperature. Since this behavior is 
not observed in the experiments of Ref.[ 7] one has to ask 
whether Eq.(16) is really applicable to the materials of 
interest. This approximation is based on the notion that 
the charge carriers move to keep each part of the host 
lattice locally electrically neutral^^, so that the Fourier 
transformed Coulomb potential, that provides the cou- 
pling between the electron and the phonon system, can 
be replaced by a constant, the deformation potential con- 
stant. However, in the systems of interest the mobile 
charge carriers are slow compared to the sound veloc- 
ity and therefore the electromagnetic potential, which 
provides the electromagnetic coupling between the elec- 
tron and the phonon system, is of very long range. Ac- 
cordingly, the electron phonon coupling constant already 
drops to zero for interaction events with very small en- 
ergy transfer. To model this effect phenomenologically 
we use the approximation 



iy{\E\) = iyoe{u; - \E\) 



(21) 



which has already been applied successfully in other non- 
equilibrium hopping problems^''. In this approximation 
the maximal amount of energy transferred in one hop is 
u). If CO is small enough, we can expand Nl{V) in Eq. 
(15) around V = V an retain olny the first term. Then 
Eq.(15) takes the simple form 



fJ-L+eU 



{jx) = evoWuj 



dVNL{V)NL{V - eU) (22) 



for eU > Lu. This equation has the same form as that 
which would be obtained for purely elastic transitions, 
although energy is exchanged with the phonon system. 
Therefore, we call this approximation the quasi-elastic 
approximation. It is Eq.(22) which has been used in the 
interpretation of the experiments of Rcf . [ 7] . 
For a DOS of the form (17) Eq.(22) yields 



(jx) = ei^oCoWeUlN^ + —^NoN^\eU\'' 

^ r(2 + 27) ^' ' J- 



(23) 



Here r(a;) is the Gamma function^^. For large ell the 
asymptotic of this expression agrees with that of the 
conventional deformation potential approximation up to 
numbers. For small eU it differs appreciably from that. 
These differences manifest themselves in particular for 
small ui. In this case the tunneling conductance ap- 
proaches the constant value 



dU 



(24) 



for Nj\eU\'^ /Nq ^ 1. The temperature dependence of 
the tunneling conductance is in this case governed by the 
temperature dependence of Nq, and thus weaker than 
that of Eq.(19). This sets the situation in the quasi- 
elastic approximation apart from that in the conventional 
deformation potential approximation and allows to de- 
cide whether the hops in an experiment are quasi-elastic 
or inelastic. If oj is small but larger than kT Eq.(24) 
crosses over to Eq.(19) if eU becomes small compared to 
CO. If a; < kT, the same dependence as in Eq.(24) is also 
observed at [/ = 0. 

The data of Ref. [ 7] are not in line with the strong 
temperature dependence of Eq.(19) (see, e.g., Fig.l of 
Ref. [ 7] Sandow et.al. (2001)). They are, however, in 
line with the Eqs.(22)-(24). Accordingly, the hops were 
quasi-elastic. 



IV. INSULATOR-LIKE CONDITIONS 

In this sections we consider the case in which the local- 
ization length is the smallest length scale in the system 
(strongly localized regime). Accordingly, the inequality 
le ^ is not satisfied. There are not enough sites 
which can be occupied by charge carriers to screen out 
the electric field on a distance of the order of the local- 
ization length. Therefore, in the insulator-like case one 
measures essentially the DOS in the space-charge region. 

For lightly doped materials far from the metal- 
nonmetal transition the space-charge region can be quite 
large. Due to this fact there is also an electric field in- 
side the sample. Therefore, the simple approximation 
eU = eEl does not hold. Instead of this relationship we 
have eU = eUL + eEl + ellR. Here Ul {Ur) is the poten- 
tial difference across the left (right) sample. The charge 
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carriers, which arc important for the tunnehng current, 
jump from the left surface of the right sample to the right 
surface of the left sample. Doing so, they have to change 
their energy by eEl. Accordingly, — i^l = &EI. In 
order to relate the difference of the chemical potentials 
to the voltage applied to the sample we express the elec- 
tric field by U. To this end we focus on the situation 
that the time for local equilibration in the right and in 
the left sample is the smallest time scale in the prob- 
lem. This implies that also the resistance of the contacts 
is large compared to the resistance of the samples, but 
small compared to the timncling resistance. In this case 
the impact of the space charge region on the tunneling 
experiment is largest. Furthermore, we assume that the 
screening of the external electric field can be described 
within the Debye-approximation. In the context of hop- 
ping transport this approximation has been discussed, 
e.g., in the Refs. [ 20] and [21]. If we use this approxi- 
mation we obtain U = Al^E+lE. Accordingly, eEl = eU, 
where U = 111/(4:1,, + 1)- Since in this case the DOS to the 
right is related to the DOS to the left by the relationship 

Nr{V, x = 1) = Nl{V -eU,x = -L^) (25) 

we obtain for the tunnel current the expression 

IXL+eij V' 

0,) = eW{l - exp(-/3eC/)) j dV J dV 

xNL{V,x = 0)NL{V'-eU,x = -L^)u{\V-V'\). (26) 

From the practical point of view the most important dif- 
ference between Eq.(13) and Eq.(26) seems to be that 
the difference between the chemical potentials is reduced, 
and therefore U is replaced by U. Due to this replace- 
ment the range of integration in Eq.(26) is getting small if 
I- This fact renders measurements of tunneling cur- 
rents more difficult. Moreover, the exponent exp(— /3ef/), 
which turned out to be negligible in the metal-like situ- 
ation, might prove to be essential in this case. Since the 
results of the experiments of Ref. [ 7] were independent 
of the tunneling distance I we conclude that in these ex- 
periments the condition I > was met. Accordingly, in 
these experiments le was at most of the order of a few 
times the average site spacing. 

In the literature the order of magnitude of the screen- 
ing length is a point discussed controversial. In Refs. [ 
16], [ 20] and [ 21] different expressions for the screening 
length have been obtained. The fact that the tunnel- 
ing current depends on the screening length le raises the 
question whether this dependence can be used in order 
to obtain further informations on Ig experimentally. We 
would like to mention that screening effects in tunneling 
experiments have been also discussed in Ref. [ 11] . 



V. CONCLUSIONS 

Starting from the usual rate equations for hopping 
transport in the impiuity band of a doped semiconduc- 
tor we have derived an expression for the tunnel current 
across the gap of a break junction device in which the 
contact material is a doped semiconductor. The fact that 
the tunnel resistance in a break-junction tunneling exper- 
iment is much larger than the resistance of the material 
leads to a separation of time scales between the tunneling 
and the dynamics inside the contact. Therefore a sepa- 
rate equilibrium inside the contacts is established with 
different chemical potentials. This simplifies the result- 
ing expression for the tunnel current as opposed to a situ- 
ation in which the contacts would be in equilibrium with 
each other and in which the tunnel and sample dynam- 
ics would be part of a common optimization or percola- 
tion problem. Due to the separation of time scales the 
situation in break-junction tunneling experiments is not 
percolation-like. The resulting expressions for the tunnel 
current look very similar to those in conventional tun- 
nel or point contact spectroscopy. They become equal 
to these expressions if metallic-like conditions apply, i. 
e. if the screening length is the smallest length scale in 
the problem. However, in the impurity band of lightly 
doped insulators the localization length is the smallest 
length scale. Therefore, the relevant contact densities of 
states are those in the space charge region. An increas- 
ing extent of the space charge region leads to a reduction 
of the difference between the local chemical potentials, 
which affects the measurement if I is smaller or of the 
order of le- If the break junction separation I is larger 
than the screening length the influence of space charge 
effects become negligible. 

We have investigated our expression for the tunnel- 
ing current in two approximations, in the conventional 
deformation potential approximation and in an approx- 
imation, which only takes into account hops with small 
energy transfer. The latter is called the quasi-elastic ap- 
proximation. In the conventional deformation potential 
approximation the tunneling conductance has a power 
like current- voltage characteristic for Pell > 1. Accord- 
ingly, the tunneling conductance scales to zero with de- 
creasing voltage in this regime. At /Self < 1 this trend 
is changed. The zero bias tunneling conductance is fi- 
nite, even if the density of states vanishes at the Fermi 
energy. Its temperature dependence is governed by the 
temperature dependence of the density of states and by 
the temperature dependence of the width of the strip of 
accessible sites. 

In the quasi-elastic approximation the expression for 
the tunneling current takes the same form as for a metal. 
For large voltages the current-voltage characteristic has 
the same asymptotic in this approximation as our expres- 
sion for the conventional deformation potential approxi- 
mation. For small voltages the quasi-elastic approxima- 
tion reflects in an ohmic tunneling conductance, which 
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only crosses over to the results of the conventional defor- 
mation potential approximation if the characteristic in- 
elastic energy is large compared to the thermal energy. In 
the opposite case it leads to a zero-bias tunneling conduc- 
tivity which depends on temperature only via the density 
of states. 

The characteristic features of the tunneling conduc- 
tance in deformation potential approximation, in partic- 
ular the strong dependence of the zero bias tunneling 
conductance, are not observed in the experiments. The 
measurements are, however, in line with our results for 
the quasi-elastic approximation. Therefore, we conclude 
that obviously only hops with very small energy transfer 
were important in the experiment. 

Let us now discuss the previous theoretical work con- 
cerning break-junction tunneling between materials in 
the hopping regime. In our opinion the conclusions^ that 
the tunneling CTirrent should be strongly fluctuating and 
strongly voltage dependent for large voltages have two 
reasons: Firstly it has been assumed that the current 
limiting hop across the tunnel gap leads upwards in en- 
ergy in contrast to our plausible reasoning. Secondly it 
was assumed that only a few tunneling events contribute 
to the current, whereas in a realistic situation the number 
N of "initial sites" for these events is very large. We have 
demonstrated that in this situation the fluctuations of the 
single current contributions do not significantly affect the 
measured current because of relation (11). Accordingly, 
we conclude, in contrast to Ref. [ 9], that statistical 
fluctuations of the tunneling current are negligible if the 
contact area of the break-junction is of macroscopic size, 
as it was, e.g., the case in Ref. [ 7]. This conclusion is 
in line with the results of the experiments of Ref. [7], in 
which the current did not show measurable fluctuations. 

Our expression for the tunneling current in deforma- 
tion potential approximation agrees, however, with that 
of Ref. [ 10] for not too small U. For very small U 
the exponent of our result differs from that of Ref. [ 
10] in two ways. First, the expression for the tunneling 
current in Ref.[ 10] yields zero for the tunneling conduc- 
tance at zero bias. This is in contrast to Eq.(19) which 
is non-zero. The reason for this difference is that in Ref.[ 
10] the occupation numbers have been replaced by step 
functions, and jumps upward have been ignored. These 
approximations become inapplicable at zero bias. They 
ignore that for kT > ell the width of the strip of pos- 
sible initial and final sites is not governed by eU but by 
kT, and that for [/ ^ upward hops with very small 
energy transfer are as likely and as frequent as down- 
ward hops. Second, the expression of Ref.[ 10] yields for 
small U a current- voltage characteristic that differs from 
our approach. The reason for the difference is that in 
our treatment we assume that charge carriers on the left 
side do not affect charge carriers on the right side. Ac- 
cordingly, the common DOS can be replaced by a simple 
product of the DOS. In Ref. [ 10], however, it has been 
assumed that the Coulomb interaction between the left 
sites and the right sites is important, and that therefore 



also the common DOS cannot be replaced by a simple 
product. We expect that such correlation effects become 
unimportant with increasing sample separation. Tunnel- 
ing experiments, however, are performed in such a way 
that the results are independent of the sample separation. 
Therefore, correlations should be not essential. 



APPENDIX 

Performing the integral (15) with the deformation- 
potential function (16) and the electronic DOS (17) we 
find that 

Ox) = eTyj.o(2/i.sa)3+27iv2A3 (27) 
x[A2/oo(A2) + 2AAtJo^(A2) + a2t/^^(A2)], 

where 

2^ r(i + a)r(i + /3) 

X3i^2(l, 3/2, 4; 2 + (a + /3)/2, 5/2 + (a + /3)/2; -A^), 

A = eU/i2hsa) and A = No/{N.,{2hsay) {.3F2 is the 
hypergeometric function). Accordingly, A is determined 
by the voltage, and A is a measure for the depth of the 
dip of the density of states. 

To get an expression for for the tunneling current for 
small A we expand Eq.(28) with respect to A. Doing so, 
we obtain 

(j,) oc CA3+« (1 - BX^ + 0(A4)). (29) 
Here C = 27, S = 12/((2 + 7)(5 -|- 27)), and 

(r(i + 7))' 



C = 



r(4 + 27) 

for A/X"/ < 1, and C = 0, B = 6/5, and 



(30) 



C = -eWiyo{2hsaf+'^'''N^A'^ (31) 



for A/X'f » 1. 

For large A we obtain 

0,) oc DX'+^^l - E/X + 0{1/X^)), (32) 

where E = 1 and 

D = ^eWM'^fisaf^^''N^A'^ (33) 
for A/X^ > 1, and = 27 -H 1 and 

(r(i + 7))' 



D = 



6r(2 + 27) 



(34) 



for A/X"' < 1. 
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In order to get some feeling for typical values of the 
parameter A we use the parameters of Ref. [ 7] . In these 
experiments voltages up to 8 mV have been used. If we 
use 2000 m/s as estimate for the sound velocity and a 
value of 2a « 10® m~^ we find that in these experiments 
the parameter A changed from 0...80. However, the data 
of Ref. [ 7] also show that in the most interesting region 
the parameter A took on only values of the order of 10 
and smaller. Accordingly, A is probably neither small nor 
large in the most interesting region in an experiment, so 
that in many cases the expression (28) has to be used for 
the interpretation of data. 
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to be taken into account since it is large compared to unity. 
To investigate the current in this case we can use the fact 
that the transition probabilities satisfy the principle of de- 
tailed balance to exchange the roles of the left and the 
right sample. Doing so, we see that the current changes 
its sign if the sign of the external electric field is changed. 
The expression for the current obtained in this way is the 
same as in Eq.(13) apart from that the sign has changed 
and from that L is replaced by R. Below we consider only 
the situation depicted in Fig.l. Accordingly, eU > in our 
calculation. 

see Ref. [ 13], p. 84, formula (4.2.18). 
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Figure 1 

Geometry of the break junction used 

in our theoretical treatment, which is 
schematically the sample geometry of 
the tunneling experiments in Ref. [ 7]. 
The hatched area is the active tunnel- 
ing region. The electric field is directed 
along the positive x-axis. In the ex- 
periments of Ref. [ 7] I/x ~ 3 — 4 



0.8 mm, L, = 1 mm and 



L = 0.2 - 0.4 mm 
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